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Compatible  Growth  and  Yield  Equations 
For  West  Virginia  Yellow-Poplar 

BRYCE  E.  SCHLAEGEL  and  D.  L.  KULOW 

Yellow-poplar  (Liriodendron  tulipifera  L.)  is  the  dominant 
species  in  its  second  growth  "yellow-poplar  type"  forest  and  in 
the  yellow-poplar  —  white  oak  —  red  oak  type  forest.  These  two 
forest  types  occupy  approximately  one-fifth  of  the  commercial 
forest  area  of  West  Virginia,  and  contain  20  per  cent  of  the  saw 
timber  volume  in  the  State  (7).'  Accurate  prediction  of  growth 
and  yield  in  these  stands  in  an  integral  part  of  effective  timber 
management. 

Growth  and  yield  are  related  variables.  The  summation  of 
the  predicted  growth  values  to  a  particular  age  should  produce 
the  yield  at  that  age.  This  has  not  been  true  previously  because 
of  the  difficulty  of  finding  the  mathematical  relationship  be- 
tween these  two  variables.  The  purpose  of  this  bulletin  was  to 
find  a  suitable  model  for  growth  that  would  give  compatible 
yield  results.  The  growth  models  were  fitted  to  the  data,  and 
tested  to  see  if  such  a  mathematical  relationship  does  exist. 


REVIEW  OF  LITERATURE 

The  first  yield  tables  for  even-aged  stands  in  the  United 
States  were  constructed  in  1894  for  white  pine  by  Gifford  Pinchot 
and  Henry  S.  Graves  (6).  These  tables  were  made  by  fitting  a 
freehand  curve  to  data  for  fully  stocked  or  "normal"  stands. 
Bruce  (1)  recognized  the  problems  involved  in  fitting  series  of 
freehand  curves  and  suggested  the  use  of  anamorphasis.-  In 
1926,  he  constructed  straight-line  yield  table  curves  of  volume 
over  age  by  site  classes  (2) .  He  found  the  curves  to  be  of  the  same 
shape  for  all  site  index  classes.  Reineke  (10)  modified  Bruce's 
method  by  using  alignment  charts  instead  of  the  anamorphasis 
curves. 

'Numbers   in   parentheses   refer  to   literature  cited. 

-Anamorphasis    is    a    freehand    method    used    to    convert    curvilinear    graphs   to    straight 

line  functions. 


D wight  (5)  suggested  the  use  of  growth  charts  in  place  of 
the  yield  table  method  to  predict  the  yield  of  even-aged  stands. 
He  used  present  diameter  as  a  basis  to  find  future  diameter  and 
volume  of  each  individual  tree.  In  the  yield  table  method,  future 
yield  is  based  upon  the  present  cubic  or  board-foot  volume  of 
the  stand.  He  recognized,  however,  that  the  errors  in  the  growth 
chart  method  may  be  large,  since  it  is  assumed  that  trees  of  a 
certain  size  at  a  given  age  will  continue  to  grow  at  approximately 
the  same  rates,  even  in  stands  of  different  densities. 

Bruce  and  Reineke  (3)  were  among  the  first  to  combine 
graphics  and  statistics  and  produced  alignment  charts  suitable 
for  yield  tables.  The  graphics  allowed  flexibility,  while  the  sta- 
tistics provided  mathematical  accuracy.  The  process  of  Bruce 
and  Reineke  involved  calculating  a  multiple  regression  equation 
and  then  converting  that  equation  into  an  alignment  chart.  This 
method  presupposes  that  the  coefficient  of  variation  of  yield  for 
the  stand  is  the  same  at  all  ages.  However,  Osborne  and  Schu- 
macher (9)  showed  for  "normal"  stands  of  red  gum  that  the 
coefficients  of  variation  for  heights  of  dominant  trees  in  the 
stand,  the  number  of  trees,  and  the  volume  are  dependent  upon 
the  age  of  the  stand.  They  demonstrated  a  method  of  construct- 
ing yield  tables  that  were  dependent  upon  the  coefficient  of 
variation  at  each  age  class.  They  found  that  the  relation  of  the 
coefficient  of  variation  of  yield  to  stand  age  can  be  used  to  de- 
termine the  form  of  the  growth  curve  of  any  site  index  from 
the  growth  curve  of  the  average  site  index. 

Until  1937  it  was  felt  that  the  "normality"  concept  was  the 
ultimate  in  forest  management.  A  few  foresters  realized  that  the 
attainment  of  a  "normal"  forest  was  impractical,  if  not  impossi- 
ble. MacKinney,  Schumacher,  and  Chaiken  (8)  used  correlation 
analysis  to  find  the  mathematical  relationship  between  the  de- 
pendent variable  volume  and  the  independent  variables  of  age, 
site,  and  stand  density.  They  produced  the  first  yield  tables  for 
"non-normal"  loblolly  pine  stands.  However,  the  summation  of 
the  growth  values  failed  to  produce  the  yield  value  given  in  the 
table  for  a  predicted  age.  In  1963,  Clutter  (4)  developed  "com- 
patible" growth  and  yield  equations.  His  work  was  based  on  the 
mathematical  dependency  that  must  exist  between  growth  and 
yield.  With  such  equations,  it  is  possible  to  predict  yield  (a  rela- 
tively difficult  variable  to  measure)  from  growth  (a  relatively 
easy  variable  to  measure) .  His  procedure  is  used  in  this  study 
on  the  first  measurements  of  yellow-poplar  growth  and  yield. 


METHODS  AND  PROCEDURES 

The  data  used  in  this  study  was  collected  in  1964  and  1965 
from  123  permanent  sample  plots  located  in  Greenbrier,  Monon- 
galia, Preston,  Randolph,  and  Nicholas  counties.  West  Virginia. 
The  plots  ranged  in  size  from  one-tenth  to  one-quarter  acre,  with 
an  isolation  strip  of  one-half  chain.  The  composition  of  the  over- 
story  consisted  of  at  least  75  per  cent  yellow-poplar  that  varied 
in  age  from  26  to  80  years  and  on  sites  of  57  to  110  feet. 

All  yellow-poplar  over  one-half  inch  and  all  other  species 
over  4.5  inches  were  tallied  by  species  and  diameter.  A  sample 
of  these  trees  on  each  plot  was  measured  for  dbh,  height  to  an 
eight-inch  top,  height  to  a  four-inch  top,  total  height,  crown 
class,  age,  bark  thickness,  and  width  of  the  last  ten  years  of 
growth.  Plot  volume  was  determined  by  the  local  volume  table 
method. 

Since  only  the  initial  measurement  has  been  made  in  this 
study,  plot  growth  was  determined  via  the  stand  table  projection 
method,  using  the  last  ten  years  radial  growth.  All  of  the  ad- 
vantages and  disadvantages  of  this  procedure  are  assumed.  A 
number  of  yield  models  were  selected  from  the  literature  (Ap- 
pendix A)  and  fitted  to  the  data.  All  models  showed  a  significant 
"fit"  to  the  data,  but  only  the  one  with  the  largest  R-squared 
value  was  used  in  this  study.  This  was  the  same  one  used  by 
Clutter  (1963).  This  yield  model  was  differentiated  with  respect 
to  age  to  give  a  growth  model  representing  volume  change  at  a 
moment  in  time.  This  model  was  then  fitted  to  the  results  of  the 
stand  projection  method  and  tested  for  "goodness  of  fit."  This 
growth  equation  was  then  integrated  into  a  prediction  equation 
for  yield.  By  summing  this  equation  over  periods  of  time,  com- 
patible yields  were  obtained.  The  steps  in  this  procedure  are 
shown  in  Appendix  B. 


RESULTS  AND  DISCUSSION 

Using  least  squares  techniques  on  model  forms  resulted  in 
a  set  of  five  equations  that  can  be  used  to  establish  growth  and 

■'The    relationship    between    height    and  diameter    for    the    sample    trees    was    found    by 

regr3ssion.    Volumes    for   selected    diameters  and    calculated    heights    were    read    from    Schnur 

(1937).    The    relationship    between    volume  in    cubic    feet    inside    bark    for    the    entire    stem 
nnd  diameter  was  then  calculated  as: 

Volume    =    4.9227   —    l.f;455D    +    0.2596D- 


yield  relationships  present  in  data  from  the  first  measurements 
of  yellow-poplar  plots.  These  equations  are: 
Cubic-foot  yield; 

LnV  =  4.7123  +  0.0071  (S)  +  0.6167  (LnB)  —  7.7335  (A' )  (2) 
Basal  area  growth  rate; 

dB/dA  =  B  (5.9358  -  LnB)  A'  (8) 

Cubic-foot  growth  rate; 

dV/dA  =  V(  1.5261  —  0.0380  LnB)A-i  (9) 

Projected  basal  area; 

LnBp  =  5.9358  —  AoAp'^  (5.9358  —  LnBo)  ( 10) 

Projected  cubic-foot  volume; 

LnVp  =  LnVo  +  1.3003  (Ln¥i-  —  Ln¥..) 
+  Ao(Ap-i  —  AoO    (0.0380  LnB..  —  0.2258)  (12) 

where: 

LnV  =  the  natural  logarithm  of  inside  bark  cubic-foot  plot 
volume, 

S  =  site  index  in  feet, 

B  =  stand  basal  area  in  square  feet, 

A  =  stand  age  in  years, 
dB/dA  =  instantaneous  rate  of  basal  area  growth, 
dV  /dA  =  instantaneous  rate  of  cubic-foot  growth, 

Bo  =  initial  basal  area, 

Ao  =:  initial  age, 

Bp  =  projected  basal  area, 

Ap  =  predicted  age, 

Vp  =  projected  cubic-foot  yield  per  acre,  and 

Vo  =  initial  cubic-foot  yield. 

Graphs  of  the  five  equations  have  been  plotted  in  Figures 
1  to  6.  Figures  1  and  2  show  yields  per  acre  for  various  age  and 
basal  area  classes  for  site  indices  of  75  to  95  (Equation  2) .  These 
are  essentially  empirical  yield  tables,  where  the  present  stand 
yield  can  be  solved  for  any  combination  of  age,  site,  and  basal 
area.  Figure  3  shows  the  basal  area  growth  rate  (Equation  8)  by 
age  and  basal  area  classes.  It  should  also  be  noted  that  the  basal 
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area  growth  per  acre  increases  with  basal  area  until  a  basal  area 
of  140  square  feet  is  obtained.  At  that  point  the  growth  rate  be- 
gins to  decline.  Figure  4  shows  the  annual  cubic-foot  growth 
rates  that  can  be  expected  at  various  ages  and  basal  areas  for 
an  initial  volume  of  3,400  cubic  feet  (Equation  9).  Although  any 
initial  volume  could  be  used  in  this  equation,  3,400  cubic  feet 
was  used  since  it  most  nearly  represents  the  present  mean  cubic- 
foot  yield. 

Residual  values  for  both  growth  and  yield  were  calculated 
and  plotted  over  age.  For  all  cases,  the  distribution  of  points 
were  random  enough  to  indicate  no  increase  in  variance  with  age. 

Solutions  of  the  basal  area  projection  equation  (Equation 
10)  appear  in  Figure  5.  These  curves  show  the  expected  basal 
area  at  various  ages  and  initial  basal  areas  by  means  of  the 
growth  rates  shown  in  Figure  3.  The  curves  for  projected  cubic- 
foot  volume  are  shown  in  Figure  6.  The  values  plotted  are  found 
by  solving  Equation  12  for  an  initial  volume  of  3,400  cubic  feet 
at  various  ages  and  basal  areas. 
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SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  study  was  to  find  compatible  growth 
and  yield  equations  for  yellow-poplar.  A  number  of  existing  yield 
models  were  fitted  to  field  data  from  123  sample  plots  and  were 
tested  for  best  "fit."  The  equation  that  accounted  for  the  largest 
percentage  of  the  variation  about  the  mean  yield  was  selected 
for  use  in  this  study.  Using  this  yield  model,  a  tentative  cubic- 
foot  growth  model  was  obtained  by  differentiation.  To  find  an 
expression  for  basal  area  growth,  a  model  for  basal  area  was  dif- 
ferentiated. Both  the  basal  area  and  cubic-foot  growth  models 
were  fitted  to  the  data  and  regression  coefficients  solved  for  by 
least  squares  techniques. 

Basal  area  and  cubic-foot  projection  formulas  were  then  ob- 
tained by  integration  of  the  basal  area  and  cubic-foot  growth 
equations.  Using  these  projection  formulas,  total  cubic-foot  yield 
may  be  solved  for  any  initial  volume,  age,  or  basal  area. 

This  study  resulted  in  the  following  equations: 

Cubic-foot  yield; 

LnV  =  4.7123  +  0.0071  (S)  +  0.6167  (LnB)  —  7.7335  (A^  (2) 
Basal  area  growth  rate; 

dB/dA  =  B  (5.9358  —  LnB)  A-^  (8) 

Cubic-foot  growth  rate; 

dV/dA  =  V(  1.5261  —  0.0380  LnB)A-i  (9) 

Projected  basal  area; 

LnB.  =:  5.9358  —  AoAi-H5.9358  —  LnBo)  (10) 

Projected  cubic-foot  volume; 

ft  ft 

LnV.'  =  LnVo  +  1.3003  (Ln?.  —  Lnfo) 

+  Ao(Ap-i— Ao-i)  (0.0380  LnB.  — 0.2258)  (12) 
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where: 

LnV  =  the  natural  logarithm  of  inside  bark  cubic-foot  plot 
volume, 
S  =  site  index  in  feet, 
B  =  stand  basal  area  in  square  feet, 
A  =  stand  age  in  years, 
dB  dA  =  instantaneous  rate  of  basal  area  growth, 
dV  dA  =  instantaneous  rate  of  cubic-foot  growth, 
Bo  =  initial  basal  area, 
A-i  =  initial  age, 
Bi'  =  projected  basal  area, 
AAi'  =  predicted  age, 
Vi'  =  projected  cubic-foot  yield  per  acre,  and 
V"  =  initial  cubic-foot  yield. 

Although  these  equations  are  mathematically  dependent,  it 
should  be  noted  that  only  one  of  the  five  significant  equations 
was  tested  (Appendix  A  and  B).  Any  of  the  others  might  have 
improved  the  compatible  efficiency  by  improving  the  growth 
model.  Then,  too,  no  attempt  was  made  to  consider  any  of  the 
other  possible  relationships  between  the  dependent  variables  and 
the  independent  variables  of  age,  site,  and  density.  This  study 
considered  only  normal  forest  and  not  such  things  as  maximum 
and  minimum  densities,  and  various  management  regimes.  These 
will  be  investigated  in  the  future  from  the  results  of  thinnings 
made  in  the  established  plots. 
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APPENDIX  A 

Selection  of  a  Yield  Model 

The  following  models  were  selected  from  the  literature  for 
use  in  this  study.  All  were  fitted  to  the  data  and  tested  by  analy- 
sis of  variance.  As  can  be  seen,  any  one  of  them  could  have  been 
used,  but  only  the  one  with  the  largest  R-  value  was  considered 
(Appendix  B). 


Yield  Model 


R^ 


Standard 
Error    Citation 


LnV  =  h^>-*  +  b,(S) 

+  h,,(A-')='* 


+  h,(LnB) 


99.9838       0.0095' 


(a) 


LnV  =  bo*=^  +  b,(A')*=^  +  b,(S)*=== 
+  b.CSA^"^ 


99.9617       0.0146'        (b) 


LnV  =  b"^==^  +  b,(S)*=^  +  b,(B)** 

+  b,(AO**  +  b,(SA-i)"^ 
+  b,(BA-^)"'' 


99.9775       0.0134'        (c) 


+  b,(BH)** 


99.6563     18.4238        (d) 


V  =  bn**  +  b,(S)**  +  b.(A)=^* 
+  b,(H)**  +b,(D)*  +  b,,(B)=^=* 
+  b,(HB  100)** 


98.7959     34.4842        (e) 


*   Significant   at  the  5  per  cent  level. 
**  Significant  at  the  1   per  cent  level, 
ns  Not  significant  at  the  5   per  cant   level. 
'  Logarithmic  scale. 

(a)  Clutter,   .Jerome  L. 

1963.  Compatible  growth  and  yield  models  for  loblolly  pine.  Forest  Sci.  9(31  :354- 
371. 

(b)  Schumacher,    F.   X. 

1939.  A  new  growth  curve  and  its  application  to  timber  yield  studies,  .lour. 
Forestry   37  :819-820. 

(c)  Olson,    David    F.,   Jr. 

1959.  A  working  plan  for  a  study  of  the  effect  of  stand  density  upon  growih 
and  yield  of  pure  yellow-poplar  in  the  southern  Appalachian  Mountains. 
Line   Project   FS-1-57-2.   Asheville   Research   Center. 

(dl     Spurr,   S.   H. 

1952.     Forest   Inventory.   The   Ronald   Press  Co.,   New  York.   476   pp. 

(e)    Smith,  .1.  H.  G.  and  Ker,  J.  W.  ^       , 

1959.  Empirical  yield  equations  for  young  forest  growth.  British  Columbia  Lum- 
berman.   Sept. 
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APPENDIX  B 

Development  of  the  Primary  Equations 

The  following  procedure  is  used  to  obtain  compatible  growth 
and  yield  equations.  This  is  the  technique  used  by  Clutter  (1963) 
in  loblolly  pine  and  is  reproduced,  in  part,  here. 

Volume  yield  —  From  the  Appendix  A,  the  most  significant  yield 
model  is: 

Ln(Vol.)  =  bo  +  bi (Site)  +  b.(Ln  Basal  Area)  +  h,(A-')      (1) 

Following  the  regression,  the  equation  became: 

LnV  =  4.7123  +  0.0071  (S)  +  0.6167  (LnB)  —  7.7335  (A^     (2) 

Noting  that  basal  area  is  a  function  of  age,  equation  (1)  may  be 
differentiated  into: 

dV/dA  =  b.VBHdB/dA)  —  b,V(A-^)  (3) 

where: 

dV/dA  =  instantaneous  growth  in  cubic  feet,  and 
dB/dA  =  instantaneous  growth  in  basal  area. 

Basal  area  growth  —  Although  Equation  3  will  provide  for  the 
cubic  foot  growth  rate,  the  term  dB/dA,  for  instantaneous  growth 
in  basal  area,  is  not  readily  available.  Therefore,  the  following 
model  was  used  to  define  basal  area: 

LnB  =  Co  +  CiS  +  CoA'^  +  C3(LnD)A-i  +  c.SA"^ 

or, 
LnB  —  Co  —  CiS  =  c,A-i  +  C3(LnD)A-^  +  c.SA"^  (4) 

where  the  symbols  are  the  same  as  in  Equation  (1)  and  D  is 
arbitrarily  the  beginning  basal  area  at  age  25.  Differentiating 
Equation  (4)  with  respect  to  age  gives: 

B'MdB/dA)  =  -AHc.oA'  +  C3(LnD)A-^  +  c.SA'^  (5) 

Substituting  Equation  (4)  into  Equation  (5)  and  rewriting  pro- 
duces: 

dB/dA  =  -B(LnB)A-i  +  CoBA"^  +  c^BSA"^  (6) 

Substituting  the  right  side  of  this  equation  into  Equation  (3) 
gives  the  cubic-foot  growth  expression: 

dV/dA  =  -b.V  (LnB)A-^  +  b.CoVA"^  +  bx.VSA"^  —  bsVA"^    (7) 
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Equation  (6)  was  fitted  to  the  basal  area  growth  data  with 
(dB/dA  +  B(LnB)A"^)  as  the  dependent  variable.  The  second  in- 
dependent term  BSA',  did  not  contribute  significantly  to  the 
regression  according  to  analysis  of  variance  so  the  quantitative 
expression  is: 

dB/dA  =  BA  '(5.9358  —  LnB)  (8) 

This  accounts  for  99.61  per  cent  of  the  variation  about  the  mean 
basal  area  growth  per  acre,  with  a  standard  error  of  0.12  square 
feet. 

Volume  growth  —  With  the  development  of  the  basal  area  growth 
model,  it  is  now  possible  to  develop  the  expression  for  cubic  foot 
growth.  Equation  (7)  was  fitted  to  the  plot  growth  data.  The 
variables  VSA"^  and  VA  '  were  non-significant  and  removed  from 
the  model.  After  analysis  this  is: 

dV/aA  =  VAM 1-5261  —  0.0380  LnB)  (9) 

and  accounts  for  93.81  per  cent  of  the  variation  about  the  mean 
cubic  foot  growth  per  acre,  with  a  standard  error  of  2.92. 

Development  of  the  Predicting  Equations 

Basal  area  growth  —  Equation  (8)  is  the  quantitative  expression 
for  instantaneous  growth  of  basal  area.  To  change  this  to  a  pre- 
diction equation,  it  is  first  solved  for  B: 

B^dB  =    (5.9358  —  LnB)A-'dA 

then  integrate  on  the  left  from  an  initial  basal  area  (B.-)  to  a 
predicted  basal  area  (Bp)  and  on  the  right  from  initial  age  (Ao) 
to  predicted  age  (Ai) . 

—Ln  (5.9358  — LnB. )  +  Ln  (5.9358  —  LnB.)  =  LnA.  —  LnAo 
and  taking  the  antilog: 

LnBi   ^  5.9358  -  A^A.  '(5.9358  —  LnB..)  (10) 

Volume  growth  —  Solving  Equation  (9)  for  V  gives: 

VdV  =  (1.5261  —  0.0380  LnB)A  'dA  (11) 

Substituting  Equation  (10)  into  Equation  (11)  (using  the  nota- 
tion B  and  A  for  Bi-  and  Ai)  produces: 

V'dV  =  1.3003  A^dA  +  0.2258  A-A  dA 

—  0.0380  AoA  MLnBo)dA 
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Integrating  on  the  left  from  initial  volume  (Vo)  to  predicted  vol- 
ume (Vp)  and  on  right  from  A^  to  Ai  (initial  to  predicted  age) 
results  in  the  following  prediction  equation  for  volume: 

LnVr  =  LnV.>  +  1.3003  (LnAp  —  LnA.) 

+  Ao(A.-^  —  Ar^)   (0.0380  LnB 0.2258)     (12) 
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